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ABSTRACT 

We present the N, O, F and Na abundance and ^^G/^^C isotopic ratio measurements or upper hmits 
for a sample of 10 C-rich, metal-poor giant stars, eight enhanced in s-process (CEMP-s) elements and 
two poor in n-capture elements (CEMP-no). The abundances are derived from IR, K-band, high- 
resolution CRIRES@VLT spectra obtained. The metallicity of our sample ranges from [Fe/II]=— 3.4 
to —1.3. F abundance could be measured only in two CEMP-s stars. With [F/Fe]=0.64, one is mildly 
F-overabundant, while the other is F-rich, at [F/Fe]=1.44. For the remaining eight objects, including 
both CEMP-no in our sample, only upper limits on F abundance could be placed. Our measurements 
and upper limits show that there is a spread in [F/C+N] ratio in CEMP-s stars as predicted by theory. 
Predictions from nucleosynthetic models for low-mass, low-metallicity Asymptotic Giant Branch stars, 
account for the derived F abundances, while the upper limits on F content derived for most of the 
stars are lower than the predicted values. The measured Na content is accounted for by AGB models 
in the 1.25 to 1.75 Mq range, confirming that the stars responsible for the peculiar abundance pattern 
observed in CEMP-s stars are low-mass, low-metallicity AGB stars, in agreement with the most 
accepted astrophysical scenario. We conclude that the mechanism of F production in current state- 
of-the-art low-metallicity low-mass AGB models needs further scrutiny and that F measurements in 
a larger number of metal-poor stars are needed to better constraint the models. 
Subject headings: stars: abundances — stars:atmospheres — stars: AGB and post-AGB — stars: 
chemically peculiar 



1. INTRODUCTION 

The understanding of the stellar nuclear production 
sites and evolution scenarios is of great current interest. 
A large number of high- resolution spectral studies target- 
ing metal-pqor obj ects selected by mining the HK survey 
(|Beers et all 11991 and the HES (IChristlieb et al.l[200l[l 
have measur ed a wide set of c hemical elements ranging 
from Li (e.g. iSpite et al.l[2005l ) to the n-capture species 
(e.g. ISneden et al.l 120031 ) . providing constraints on how 
the stars formed as well as on the nucleosynthetic pro- 
cesses that took place in the early Galaxy. 

In contrast, only one single measurement of F in ex- 
treme ly metal-poor (EMP) stars has been published to 
date (jSchuler et al.l 120071 ). The most accessible lines 
suited for F measurements, the vibration-rotation tran- 
sition of the HE molecule, are located in the K band, 
requiring the use of high-resolution, IR spectrographs 
mounted on 8 m class telescope to target halo stars, a 
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configuration not available until recently. 

Fluorine is an element of particular interest, extremely 
sensitive to the physical conditions within stars. Al- 
though the single stable F isotope, ^^F, is not involved 
in the main reactions taking place in the cores of stars, 
it can be created and destroyed in several different ways, 
and its dominant source is not yet clear. Theoretical 
modeling has indicated (at least) three possible sites for 
F production. 

iWooslev fc HaxtonI (|1988| ) proposed first that F is pro- 
duced by neutrino spallation on ^°Ne; ^°Ne(i^,i/'p)^^F. 
Second, modeling shows that low- metallicity, low- 
mass (M<3-4M0) AGB st ars synthesize very large 
quantities of F (see e.g. iCristallo et al.l l2009( ) via 
the reactions i'^N(n,p)"C(a,7)i80(p, a)i5N(a, 7)19F 
and i4N(a,7)i8F(/3+)i8o(p,a)i5N(a,7)i9F, where 
the neutrons are provided by ^'^C(a, n)^^0 and the pro- 
tons mainly by ^'*N(n,p)^^C (see iJorissen et aflll992|) . 
On the other hand, massive AGB stars are expected 
to destroy F via Hot Bottom Burning ( HBB) (see 
e.g. iSmith et alT 2005; Karakas fc Lattanziol 12007). Fi- 
nallv. jMevnet &i Arnould 2Q0Qi) proposed that Wolf- 
Rayet stars are a possible source for F, through the 
i'*N(a,7)i« F(/3+)i80(p,a)i5N(a,7)i9F chain. While 
the physical conditions for F production exist in several 
phases of stellar evolution, the precise contribution of 
each of these three sites to the Galactic F evolution is 
not known. 

The measurements of F content in any object have 
been few and give con flicting evidence ab out the origin 
of F in the Universe. iRenda et al.l (|2004l ) showed that 
the inclusion of all of sites in Galactic chemical evolu- 
tion models is necessary to reproduce the F abundances 
measured in Milky Way field stars. In their model, AGB 
stars are responsible for significant amounts of F pro- 
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duction in the early Universe, because of the metalhcity 
dependence of their yields, while WR stars are signif- 
icant contribu t ors a t solar and super-solar metallicities. 
iJorissen et al.l (|1992[ ) derived F abundances in giants that 
recorded the chemical evolution (K and M giants), in 
AGB giants that were dredging up freshly minted F from 
their interiors (MS, S, SC, N and J giants), and in gi- 
ants that had been polluted by one of the former (Ba 
giants) . They suggested that the overproduction of F in 
the AGB stars (factors of 3 to 30 with respect to the 
solar system in the most extreme S and N stars) showed 
that AGB stars are likely to major contributors to the 
Galactic F abundance, at least at metallicities close to 
solar. In addition, the largest F overabundances mea- 
sured could not be explained with standard AGB mod- 
els and required additional mixing to achieve the d esired 
am ounts. However, r ecent results bv lAbia et al.l ()2009l ) 
and lAbia et al.l (|2010( ) suggest that because of a possible 
lack of proper accounting for C-bearing molecules (i.e. 
CH, CN, CO and C2) contribu tion, the F abundances 
reported in IJorissen et al.l (|1992D for solar metalhcity gi- 
ants, which relied on the same HF line used in the present 
analysis, had been overestimated. Therefore, the long 
standing problem of the discrepancy between the high 
[F/Fe]3 and [F/0] ratios measured in solar metalhcity gi- 
ants and the low ones resulting from the models, is likely 
solved by the adoption of a more complete molecul ar line 
list fo r the synthetic spectra. In fact, the Cristallo et al.l 
(j200il) models account fully for the F abundances mea- 
sured in several C-stars ranging from solar metalhcity 
to slightly m etal-poor ([Fc/H]-0.0 to —0.5) (jAbia et al.l 
I2OO9II2OIOI) . 

While the above results make it clear that solar metal- 
hcity AGB stars make at least some F, there is not yet 
a consensus on the overall contribution of AGB stars 
throughout the che mical evolution of stellar systems. 
ICunha et al.l ()2003D measured F in red giants in the 
LMC and uj Gen. They found that F/0 declines as the 
O abundance decreases and the two w Gen giants have 
particularly low F/0 values. They argue that their re- 
sults are qualitatively consistent with most F produc- 
tion coming from either neutrino nucleosynthesis or WR 
stars rather than AGB stars in these systems, which cer- 
tainly have a different chemical evolution than the Milky 
Way. In particular, the low F abundances in metal-poor 
([Fe/H]'^ —1) a; Gen stars that are enriched in s-process 
elements is difficult to reconcile with the idea that metal- 
poor AGB stars are major contributors to F production. 
F has also be en measured in stars in the Galactic bulge 
(jCunha et al.| [2"008). The F/0 resuhs there can be ex- 
plained by contributions by both AGB and WR stars, al- 
though the lack of an s-process enhancement in the most 
F-rich bulge stars suggests that AGB stars may play a 
less prominent role in the bulge than has been inferred 
for the disk. 

Measurements of flu orine in the interstellar medium 
(|Federman et al.l I2005D show no evidence of F over- 
abundances due to the neutrino process in SNII. The F 
pr oduction in rotating WR stars has been reconsidered 
bv lPalacios et al.l (|2005l ) , who found that F yi e lds ar e sig- 
nificantly lower than the iMevnet fc Arnouldl ()2000D pre- 



dictions, indicating that their contribution to the Galac- 
tic F budget would be negligible. These results suggest 
that at low-metallicity AGB stars play a major role in 
F production, although this idea would leave the w Gen 
results unexplained. The idea of AGB stars as producers 
of F is supporte d by the large F enh ancements found in 
post- AGB st ars (IWerner et al.ll2005D and planetary neb- 
ulae (see e.g. lOtsuka et al.ll2008l and references therein), 
the progeny of AGB stars as well as the results for abun- 
dances in Milky Way AGB stars cited above. However, 
these observations are focused on [Fe/H]^0 metalhcity 
stars and do not address directly the question of the pro- 
duction of F production in metal-poor AGB stars. 

While Na measurements have been obtained for a much 
larger number of metal-poor stars, no quantitative study 
has been performed so far on the role of AGB stars in 
its production, in particular at low metalhcity. In ther- 
mally pulsing AGB stars, Na is mostly synthesized by 
proton captures on ^^Ne. Models indicate that the main 
source of Na in low-mass AGB stars, at least at solar 
metalhcity, is the creation of a Na pocket, located at the 
top (or near the top) of the ^''N pocket. At lower metal- 
licities, other mechanisms, such as the neutron capture 
on ^^Ne, which can occur both during the radiative ^'^G 
burning and during the convective ^^Ne burning, may 
become important (Gristallo et al 2006b and references 
ther ein). In general, acc ording to e.g., Gristallo et al. 
and lBisterzo et all (|2010( ). the production of Na increase 
by about 1 dex or more going from solar to [Fe/H] -2.3, 
a similar metalhcity to that of the stars in our sample. 
This is due to a higher efficiency of all the Na nucle- 
osynthesis channels described above. The testing of such 
predictions is important to constraint details of the key 
reactions involved. 

Carbon-enhanced metal-poor (GEMP) stars provide 
an opportunity to directly measure the F and Na pro- 
duction in low-mass, metal-poor AGB stars. These 
stars are chemically peculiar objects, characterized by 
an overabundance of G 
20% of stars below [Fe/H 



G/Fe]>lE3) acc ounting for 10- 
< - 2.5 (M arstel ler et all [20051 : 



^ Hereafter:log e(A) 
log(Y/Y0). 



log N(A)-hl2. [X/Y] = log(X/Xo) 



iGohen et all 120051 : iLucatello et"ar.l i2006'). Less than a 
third of GEMP stars exhibit no enhancement in heavy 
elem ents (GEMP-no) , while most of these objects (over 
70% lAoki et al]|2007[ ) are characterized, by an overabun- 
d ance of n-captur e , s-pr ocess, elements (CEMP-s). 

ILucatello et al.l (|2005[ ) showed that likely all C-rich, 
extremely metal-poor stars with s-process enhancement 
(the GEMP-s stars) belong to a binary system. CEMP- 
s are then the metal-poor analog to the classical CH 
and Ba stars: low-mass stars (M~0.8Mq) whose slightly 
more massive (between ~1.2 and ~2.5Mq, the exact 
range depending on metalhcity) companion, now a faint 
white dwarf, dumped material processed during its AGB 
phase on their surfaces, leaving its chemical fingerprints 
in the composition of their envelopes. Therefore, the nu- 
cleosynthetic processes taking place in extremely metal- 
poor, low-mass (~1.5Mq) stars, now long extinct, can 
be investi gate d through the study of GEMP-s stars char- 
acteristicsPl. 

Some authors use different cutoffs values depending on the 
evolutionary state of the star, with cutoff enhancement for giants 
as low as [C/Fe]~0.5 fsee lAoki eraI|[2007D 

Because of its high Eu abundance, not fully accountable with 
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The origin of C EMP-no object s is, on the other hand, 
still a mystery. iFuiimoto et alj ()200Cl[ ) suggested that 
they may have formed as chemically "normal" low-mass 
stars and became C-enhanced through a path of self- 
enrichment due to anomalous mixing processes specific 
to low-metallicity stars. Alternatively, as iRvan et al.l 
(jlOOSi) suggested, they could have been born from C- 
rich gas, possibly polluted by a previous generation of 
supernovae whose fall-back avoids the ejection of heav- 
ier elements during the explosi on (e.g., from h igh- mass, 
rapidly rotating stars see Mev net et al. I [20061 and ref- 
erences therein). Objects like the recently discovered 
extremely-metal po or, C-rich Damped Lyman-a system 
(jCooke et al.ll20ldl) might turn out to be the connection 
between the yields of the Pop III stars and their later 
incorporation into CEMP-no stars. 

On the other hand, their abundance patterns could 
arise from early AGB transfers from low-mass stars be- 
fore a ny considera ble s- process element produc tion took 
place ()Rvan et al.ll2005l:lMasseron et al.ll2010aD . or alter- 
natively from an AGB stars whose evolution was trun- 
cated by binary interaction (see, e.g., Izzard & Tout 
2003) with its companion, the presently observed object. 

The measurement of a value or a strong upper limit on 
F abundance is of crucial importance to probe the ori- 
gin of the chemical pattern observed in CEMP-no stars 
given that F can be synthesized before the bulk of th e 
s-process-element production. In iCristallo et al] ()2009( ). 
for instance, as early as the third-dredge-up episode for 
a 2M0 at [Fe/H]=-2.3, [F/C]=-0.4 and [F/Fe]=1.65, 
[Ba/C]=-1.5and [Ba/Fe]=0.5 (and [C/Fe]=2.05). There- 
fore, a measurement in CEMP-no stars of an amount of 
F comparable to C would strongly argue in favor of AGB 
enrichme nt. 

Schulcr et al.l (|2007f ) measured the F abundance in a 
extremely metal-poor star for the first time, deriving 
an abundance of loge(F)= -1-4.96 ± 0.21 corresponding 
to an abundance ra tio |"F/Fe] = 4-2.9 fo r the CEMP-s 
star HE 1305-1-0132. 'Lugar o et all (|2008[ ) compared this 
value to existing nucleosynthesis and mass transfer mod- 
els. Conclusion that an object with such an extreme F 
content should be exceedingly rare, while most CEMP-s 
stars are expected to exhibit a noticeable but smaller F 
overabundance . 

We here present F abundance measurements for two 
CEMP-s and upper limits for eight more CEMP stars 
(six CEMP-s and two CEMP-no) and discuss the im- 
plications of our result on our current understanding of 
AGB nucleosynthesis. 

2. OBSERVATIONS 

Data were collected using CRIRE^HI (CR yogenic high- 
resolu tion InfraRed Echelle Spectrograph. iKaeufi et afl 
I2004D at VLT UTl (Antu) in a series of service runs 
between October 21th 2007 and February 28th 2008. 
The slit-width was set to 0'.'4 yielding a resolution of 

standard s-process nucleosynthesis, one of the stars in our sample, 
HD187861, has been included in the CEMP-rs category. However, 
since all the likely formation scenarios invoked to explain the abun- 
dance patterns of t hese ob jects include mass transfer from a n AGB 
companion (see e.g. lJonsell et al.l2006l : lLucateIIo et al.l2009l ). as far 
as the present discussion is concerned HD187861 can be considered 
as an CEMP-s 

^2 www.eso.org/instruments / crires / 



R = A/AA ~ 50,000, where AA refers to the width 
of the resolution elements, with a binning of 1x1. We 
used the 24/-i/i setup, which leads a wavelength coverage 
from 2299.3 to 2311.8 nm, 2315.2 to 2327.3 nm, 2330.4 
to 2342.0 nm, and 2345.0 to 2356.0 nm for chips #1 to 
4 respectively. Our sample also included a well-studied 
C-normal metal-poor star (HD 122563) to be used as a 
comparison star. The typical signal-to-noise ratio for the 
program stars at ~2335 nm is about 120 per pixel. An 
early-type, fast-rotating star was observed right before 
or right after each one of the program star observations. 
To improve the accuracy of the telluric line subtraction 
procedure, those objects were chosen so that the science 
and calibration stars were observed at similar air-masses. 
The typical signal-to-noise ratio for the telluric standard 
stars at ~2335nm is about 250 per pixel. 

The raw frames were processed with the CRIRES 
GASGANO pipeline (version 1.6.0). The ID science 
and standard star spectra were wavelength-calibrated 
separately using the numerous telluric absorption lines 
present on all of the four detector arrays. Such pro- 
cess was performed separately for the science and tel- 
luric standard star spectrum because of the limited re- 
producibility of the Echelle grating position. The science 
spectra were then divided by the appropriate standard 
star spectra to correct for the telluric lines and the illu- 
mination pattern in the best possible way. It is worth 
noting that while the telluric lines are strong enough to 
be used for wavelength calibration, they are weak enough 
to be accounted for by standard star division; hence their 
effect on the abundance measurements presented in this 
paper is negligible. 

More details about the observations can be found in 
Table m 

3. ATMOSPHERIC PARAMETERS AND ATMOSPHERIC 
MODELS 

The determination of the adopted atmospheric param- 
eters can be obtained photometrically, via the derivation 
of the effective temperature Teff from accurate Tcff-color 
relations and that of surface gravity, log g from a suitable 
isochrone. Alternatively, T^s and log 5 can be derived 
spectroscopically, through the requirement of, respec- 
tively, of excitation and ionization equilibrium. In both 
cases, the micro turbulent velocity ^ is determined from 
minimizing the trend of derived abundance and equiva- 
lent widths of the Fe lines. 

The spectral range covered by the program spectra do 
not include any Fe lines and in general not enough fea- 
tures to allow the derivation of the atmospheric param- 
eters. Therefore, we have relied on the atmospheric pa- 
rameters published in the recent literature or derived on 
the basis of photometry and our own proprietary opti- 
cal spectra (for details see Masseron et al. 2010a). The 
adopted atmospheric parameters and C and Ba abun- 
dances (which cannot be derived from our spectra as no 
suitable features are present in the K band), along with 
their uncertainties and sources are listed in Table[TJ Note 
that in some cases the uncertainties are not listed as the 
source paper did not report them. 

The chemical composition of the stellar atmosphere 
can affect its structure: the higher the content of met- 
als the larger the number of free electrons, available for 
the formation of the H" ion, the main source of contin- 
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uum opacity among stars with spectral types later than 
A. Solar-scaled model atmospheres, which are generally 
(appropriately) adopted for the analysis of stars, likely 
provide different temperature profiles from those that are 
expected to be found in CEMP stars, due to their large 
C and N enhancements. In our abundance analysis, we 
adopt C and N enhanced n iodel atmospheres belon ging 
to the OSMARCS family (iGustafsson et all [20081 ) cal- 
culated using the appropriate atmospheric parameters 
and C and N enhancements for each star. Such models 
have a steeper temperature profile than their solar-scaled 
counterparts, affecting especially the computation of the 
strong lines, which form in the outer layers. The discus- 
sion of the calculation of the C and N-enhanced mod- 
els, their features and their impact on the abundance 
determinations is beyond the aim of th e present paper. 
For mo re details the reader is referred to IMasseron et al.l 

(pnTohl) . 

4. ABUNDANCE ANALYSIS AND COMPARISON WITH 
THE LITERATURE. 

Elemental abundances for N, O, F, Na and C iso- 
topic ratios were obtained by synthetic spectra fits to 
the ID extracted spectra for the sample stars. We use d 
the TURBOSPEC synthesis code (I Alvarez PledllQQSh . 
which shares routines and input data with OSMARCS. 
The code is characterized by chemical equilibrium includ- 
ing 92 atomic and over 500 molecular species, and Van 
der Waals colli sional bro adening by H , He, a nd H2 follow- 
ing [Barkle m fc Aspelun d-JohanssonI (120051 ). Smooth- 
ing for the spectra was determining by simulta- 
neous fitting of a number of weak unblended CO 
lines. 

A partial analysis of this dataset has been pre- 
sented in Lucatello et al. (2009). The present pa- 
per presents a reanalysis of the data, providing 
measurements for a wider set of elements. It is 
based on improved reduction of the spectra (us- 
ing a different pipeline version) and makes use 
of more accurate atmospheric parameters deter- 
minations. Hence the present results supersede 
those presented in Lucatello et al. (2009). 

O was measured from the fit of several ^^C^^O lines be- 
longing to the V = 2-0 and 3-1 overtone bands, as well as 
the V = 4-2 and ^^CO i' — 2-0 band-head regions. Given 
the high C content of the stars under consideration, the 
CO lines are very sensitive to O content, hence allowing 
for a very accurate O abundance measurement. ^^C/^'^C 
isotopic ratios were measured from the simultaneous fit- 
ting of several ^^C^^O and ^^C^^O lines. O abundance 
was then re-measured with the appropriate C isotopic 
ratio and the procedure was repeated to convergence. N 
abundance measurements are obtained from a number of 
CN lines scattered throughout our spectral range. 

Three Na lines are included within the available spec- 
tral region: 2334.842 nm, 2337.895 nm and 2337.914 nm. 
The latter is the strongest and, unlike the other two, is 
not severely blended with CO features, making it the 
best candidate for Na abundance measurements. It falls 
in a range of strong telluric contamination, therefore, in 
spite of careful and accurate modeling and subtraction of 
telluric features, in four out of eleven (10 CEMP and one 
EMP) stars only upper limits (Icr) to the Na abundances 
could be placed. 



Finally, F abundances (and Icr upper limits) were 
measured from the fitting of the unblended HF line at 
2335.833 nm. Two more HF lines are present in our spec- 
tral range, 2324.062 and 2348.805 nm; however both of 
them are blended with CO features, hampering the accu- 
racy of F abundance determinations obtained from their 
fitting. Figure [T] shows the fit of synthetic spectra to the 
HF feature at 2335.833 nm. 

The line list adopted are from iGoorvitclTI (I1994D for 
CO, Kurucz (1995) for CN and Na. and IDecin I (|2000[ ) 
for HF. They were checked for consistency by fitting 
the Arcturus spectrum and that for our comparison star 
HD 122563. 

In a few cases (CS 29502-92, CS 30314-67 and 
HD 122563) no O abundance in conjunction with the lit- 
erature C abundance could yield satisfactory fits of the 
CO features. Hence we sligthly adjusted the C abun- 
dance to improve the fit of such features. This discrep- 
ancy could be due to the adoption of different model 
atmospheres in the present analysis and in the literature 
studies in which the adopted C-abundances are reported. 
In particular the inclusion of the appropriate opacities 
for the C and N enhanced atmospheres can result into 
the derivation of slightly different abundances from those 
obtained with solar-scaled models, which are generally 
adopted in the literature. 

The uncertainties in our measurements come from two 
sources (excluding systematic NLTE and 3D atmosphere 
effects, and atomic and molecular parameters uncertain- 
ties): error in the fitting of the spectral features, and due 
to the atmospheric parameters. 

In the case of the fitting uncertainty, we estimate this 
component of the error by how well fitted we can match 
the synthesis to the actual spectrum (we adopt ~0.1 dex) 
To estimate the errors due to the uncertainty in the at- 
mospheric parameters for each element, we determined 
individually the sensitivities of the derived abundances 
to the adopted Tg//, surface gravity logg, C enhance- 
ment, and [Fe/H], and then summing in quadrature the 
resulting uncertainties associated with each parameteiF^. 
Note that for the stars for which no estimates for the at- 
mospheric parameter uncertainties were provided in the 
source paper, we adopted 100 K for Teff, 0.3 in log 5, 
0.15 in [Fe/H] and 0.2 in [C/Fe]. The sensitivities to the 
adopted atmospheric parameters are provided in Table [3] 
for HE 1152-0355 and HD5223, which beside being the 
only two stars for which F is measured, are representative 
of the present sample for atmospheric parameters. The 
estimates of the errors reported are for abundances (e.g. 
log e(F) or [F/H]), however the abundance ratios between 
two species (e.g. [F/Fe]) have typically smaller errors 
since the effects due to uncertainties on the atmospheric 
parameters partially cancel out. As above mentioned, 
we cannot obtain the information about Fe sensitivity 

^•^ Strictly speaking this is only appropriate in case of indepen- 
dent errors, and in general one needs to take into account the co- 
variances. Given that the atmospheric parameters for most of the 
objects in our sample are taken from literature sources, the deter- 
mination of the covariance terms is essentially impossible. On the 
basis of the data for the three stars for which atmospheric parame- 
ters were derived from our proprietary high-resolution spectra, we 
obtained estimate the variations in log 3 and [Fe/H] for an increase 
of 200 K in Te// , finding an average of 0.5 dex and 0.15 dex respec- 
tively. Hence, T^f f is clearly the dominant source of uncertainty 
(see Table [Sj, and the covariance terms can be neglected. 
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to atmospheric parameters in our sample stars, therefore 
we cannot actually calculate the errors on the abundance 
ratios. In the following plots and discussion we adopt for 
the ratios of F, C, N, O and Na to Fe the same errors, 
keeping in mind that they are overestimates of the actual 
ones. 

The resulting N, O, F and Na abundances, along with 
the observational errors, are listed in Table H] 

The HF lines in the spectrum of the comparison, C- 
normal star HD 122563 were undetectable as expected, 
thus only an upper limit could be placed to the F abun- 
dance. This upper limit is consistent with the galactic 
chemical evolution model for such element, which pre- 
dicts a value of [F/0]~ -0.6 dex at loge (O)=7 (the 
minimum metallicity taken into account, see lRenda et al.l 
[2004), while we find [F/0]< -0.1 dex at loge(0)=6.8. 

For the program stars, as shown in Table [3J F measure- 
ments could be derived only for two out of ten program 
stars; with metallicities of [Fe/H]=— 1.2 and —2.0 they 
are the object with the highest Fe content of the sample. 
For the other eight, more metal-poor objects, only upper 
limits could be placed. Because of the high sensitivity of 
the HF molecule to temperature, as shown in Fig. [21 the 
level of upper limits on the F abundance is essentially set 
by the Tcff in the atmosphere of the observed stars. 

As F has been measured only in one (not part 
of our sam ple) extremely metal-poor star before (see 
iSchuler et al. 2007), our results for F abundances can- 
not be compared to literature ones. However measure- 
ments for N, O and Na have been previously reported 
for a subset of our sample. Fig. [3] shows the excellent 
agreement with literature results for these species, which 
are all measured from transitions in the visible spectra. 

5. DISCUSSION 

With measurements or upper limits for eight stars, we 
find a range in [F/Fe] and [F/C-l-N] ratios in CEMP stars 
(see e.g. Fig. [5]), with CEMP-s [F/Fe] values ranging 
from -1-1.44 dex to less than -1-0.6 dex. This confirms both 
that AGB stars are characterized by range of efficiency 
in its production (likely due to different metallicity and 
and mass), as predicted by theoretical low-mass, low- 
metallicity nucleosynthetic AGB stars models. The two 
stars with measured F, HD 5223 and HE 1152-0355, are 
both CEMP-s stars, and the presence of both F and s- 
process elements in their atmospheres is consistent with 
a picture where they were polluted by a low-mass AGB 
star. 

There is no object in our sample with a F con t ent as 
high (or close) to that presented in ISchuler et al.l (|2007[ ) 
for HE 1305-f0132. In fact, even the highest value for 
an upper limit in the present sample, [F/Fe]<2.5 is al- 
most a factor of three lower than those measured in HE 
1305-1-0132. Possible sources of systematic offsets include 
difference in the model atmospheres (we use spherical 
C-enhanced models, whereas Schuler et al. adopt C- 
enhanced, plane parallel models) and in the molecular 
transition data. 

In our analysis we used spherical C-enhanced mod- 
els. These are cooler at the surface than the C-enhanced 
plane parallel models adopted by Schuler et al, as can be 
seen in Figure [U which shows the difference of the ther- 
mal structure of the two models for one of our star in our 
sample (HE1152-0355). The discrepancy in T reaches 



lOOK. Abundances derived from lines originating in the 
outer layers of the stellar atmosphere are thus lower in 
the spherical model case than in the plane parallel one. 
However, such affect is expected to affect strong lines (in 
our case mostly CO lines), while its effect on weak lines 
such as HF is expected to be very small. 

As discussed in Section [3 for the HF molecular tran- 
sitions we use the data inlDecinl (l2000fl . while Schuler 
et al. use those in lJorissen et al.l ()1992f) . While the g/'s 
for the HF(l-O) R9 line (the only present in the Schuler 
et al. spectral range) and dissociation potential adopted 
are virtually identic a l, the excitation potential for the 
transition in iDecin I (|2000[ ) is 0.227eV, weU below the 
0.48 adopted by Schuler et al. Such difference leads to 
an offset in the F abundances measured from the HF(1- 
0) R9 line. To asses such difference, we derived the F 
abundances adopting the same molecular transition pa- 
rameters as Schuler et al., and found that the values ob- 
tained are at most 0.3 dex larger using the Jorissen et al.l 
()1992|) excitation potential than with the Decin I ()2000f ) 
one. Hence, even taking into account the offset due to the 
difference in adopted excitation potential, the F abun- 
dance and [F/Fe] ratio measured in HE 1305-1-0132 are 
higher than any of the measurements or upper limits de- 
rived in our sample. 

It is noteworthy t hat the Fe abun dance in 
HE 1305-f 0132 adopted in lSchuler et~all ()2007[ ) was quite 
uncertain. In fact, given the limited spectral coverage 
in the Phoenix data and the lack of Fe features in that 
range, the authors relied on [Fe/H] m easured from inter - 
mediate resolution spectra analysis ([Beers et al.ll2007t) , 
which reported an [Fe/H]=— 2.5±0.5dex. More recent, 
high-resolution optical data suggest that HE 13054-0132 
is in fact more me tal rich, with [Fe/H]= —1.92 dex 
([Schuler et al.l l2008| ) . When taking into account this 
Fe abundance the [F/Fe] ratio lowers by '^0.6dex, to a 
similar level to those measured in our sample. Moreover, 
the F abundance itself, loge(F) can also be affected as 
the adoption of an appropriate, more metal-rich model 
atmosphere, whose outer layers are cooler than those of 
the metal-poor model atmosphere adopted previously, is 
expected to result in a lower F abundance. 

Although the present sample shows that there is a 
range in [F/Fe] abundances in CEMP stars, it is not 
yet clear what the maximum F production is, and larger 
number of F abundance measurements in CEMP (and 
EMP) stars is crucial to probe the matter. 

5.1. Theoretical AGB predictions for F production 

Since iJorissen et al.l ([1992D . considerable effort has 
been expended to model its production and explore the 
effect of uncertainties in the adopted parameters. The 
main issue to address was that the nucleosynthetic AGB 
models failed to acc ount for the h igh [F/Fe] and [F/0] 
rations observed in IJorissen et al.l |l992) in metal-rich 
in trinsic and extrin sic C -stars. Whi l e the recent results 
of lAbia et ahl ([2009f ) and lAbia et"aL[ ([20100 seem to have 
reconciled observations and model predictions at solar 
and slightly subsolar metallicity, large uncertainties in 
the ^^F production in AGB stars, especially in the low- 
metallicity regim e , still exist. 

IJorissen et al.l ([1992D proposed as a main mecha- 
nism to produce F in AGB stars the nucleosynthesis 
chain ^'^N{a,jy^F{P+)^^0{p,a)'^^'Nia,j)^^F during the 
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thermal-pulse (TP) in the He intershell. In this case, 
protons for the ^^0(p,a)^^N reaction are provided by 
"'^'*N(n,p)^^C, and neutrons are provided by ^^C(a,n)-'^^0, 
where ^'^C is prese nt in the H-burning shel l ashes at 
the TP activation. iGorielv fc Mowlavl ( 2000f ) proposed 
that additional F could be made during the interpulse 
phase, where the reaction ^*N(n,p)^^C(a, 7)-^*0(p,a)^^N 
occurs, the neutrons being provided by the ^'^C pocket 
and the protons recycled from the -'^"'N(n,p)-'^"'C. The pro- 
duction of ^^F is then completed in t he next pulse by a- 
capture on ^^N. iLugaro et al.l (|2004D described the con- 
sequence of the reaction rates uncertainties and of the 
inclusion on a ^'^C pocket on the -'^^F yield, showing that 
their combined effect can affect i t by as much a s 40% 
in low-mass stars (M< 3Mq). iKarakas etldl (|2OO80 
instead focused on the uncertainties on ^^F produc- 
tion stemming from the "'^^F(a, p)^"'^Ne reaction, whose 
rate was determined experimentally only in 2006 (see 
ILee et alJ [2006) . Such a reaction was neglected in the 
previous work (|Lugaro et al.l [200^ . The effect of the re- 
action uncertainties on F production is found to increase 
with decreasing metallicity, reaching a factor of two for 
a 2M0 star at [Fe/H]= — 2.3dex. Another important 
source of uncertainty is mass loss. Stancliffe & JcfFcrv 
([2007') explored the effects of different prescriptions for 
mass loss at the beginning of the thermally pulsing- 
asymptotic AGB phase on F production, showing that 
it could change it by as much as a factor of ~5 (see also 
iStancliffe fc Glebbeek 1 12008| ) . In particular, Stancliffe 
(|2009'), without any ^'^C pocket, using a given mass 
loss pres cription, a di fferent thir d dredge-up efhcienc; 
from the ILugaro et all (|2004f ) and IKarakas etall ()200. 
models, at M=2M0 and [Fe/II]=-2.3 predicts an F 
abundance which is ~10 times lower than than the 
ILugaro et all (|200l . 

AH these three family of models mentioned, however, 
do not take into account the effects induced by carbon 
enhancement on the opacities of the cool external lay- 
ers of AGB stars. When the star becomes C-rich and 
the number or C ions exceed that of O ions, the opac- 
ity of the external layers is affected (increased) by the 
massive presence of C-bearing molecules (mostly CH, 
C2, CN and CO). This is expected to causing the en- 
velope to expand and t he star to become la rger and 
cooler (seelMarieo II2002D . iCristallo etall (I2009D (see also 
|http://www. ugr.es/~cristallo/data-online. html) calcu- 
lated models of low metallicity, low-mass AGB stars in- 
cluding the effects of C- and N-enhancements on stellar 
opacities. This has the effect of decreasing the ^^F pro- 
duction, by roughly a factor of 5 with respect to the 
standard models (i.e. obtained usin g the r ecommended 
reaction rates) from IKarakas et al.l (2008'); it is note- 
worthy that using a solar-scaled composition Cristallo 
finds the same results as the standard K arakas model. 
It sho u ld be noted that wh ile th e in the ILugaro et al.l 
(|200l . IKarakas et all (pOOl and ICristallo et al.1 (|2009[ ) 
models ^^F prod uction peaks at around 2.2 Mq, in the 
IStancliffe I (|2009D it is reached at 1.5Mq. 

5.2. Comparison with models: C, N, and C isotopic 

ratios 

In this section, the results from the previously dis- 
cussed four different families of models (which pro- 



vide predictions of mass production of various elements 
among which C, N, O, Na and F) are compared to the 
data. Under the assumption that the CEMP-s stars ob- 
served have been polluted by mass transfer from an AGB 
companion and that the amount of C, N, Na, and F in the 
original (pre mass-transfer) atmospheres are negligible 
with respect to that accreted, the [N/C-^NQ [F/C-I-N] 
and [Na/C-|-N] ratios observed in CEMP stars should re- 
flect that produced by the AGB star companions, inde- 
pendently from the amount of dilution (ratio between the 
mass accreted and that of the envelope of the observed 
star) undergone by the accreted material. Figures [5] 
and[6]show the behaviour of N and ^^C/^^C in the 
sample stars. CEMP-s stars are plotted in black 
and CEMP-no stars in green. Given the lack of 
understanding of the nucleosynthetic process to 
which the latter owe their peculiar composition, 
the inclusion in the subsequent plots of the two 
CEMP-no in the comparison to the AGB models 
is, strictly speaking, not appropriate. However, 
since one of the possible scenario is indeed anal- 
ogous to that invoked for CEMP-s stars, we in- 
clude them in the plots for completeness (using a 
different color for easy differentiation). In Fig. [5|, 
for each individual model different degrees of di- 
lution correspond to different positions along its 
correspondent line. 

The choice of the sum of the two species C-|-N 
rather than a single one at a time is because it 
is a much more solid prediction in stellar evolu- 
tion models. It is well known that the standard 
AGB models fail to reproduce the N abund ances 
obser ved in C-rich stars (see e.g. Jorissen et al.l 
119920 . We find similar results in our sample. The 
common explanation of such apparent disagree- 
ment is that extramixing processes take place at 
the bottom of the AGB envelope, where part of 
the envelope material undergoes the CN cycle, 
turning carbon into nitrogen (see the discussion 
at the end of this section). 

As can be noticed from Fig. [5| the low mass 
models predict a N abundance which is lower 
than the measured ones, which are closer to 
the predictions of more massive AGB stars (2.5- 
4M0) models, where HBB takes place. 

A similar conclusion can be drawn from Fig. [6|, 
which plots 12C/13C with respect to C-|-N. It is 
worth noting that, unlike the previous case with 
N and those of F and Na that will be discussed 
later, the initial composition of the star (before 
mass transfer) is not negligible when plotting the 
model prediction. In fact, while essentially all 
C (as ^^C) comes from mass transfer, accord- 
ing to the low-mass AGB models (which do not 
include extramixing processes) the ^^C/^'^C ra- 
tios are very high, thus according to these mod- 
els, very little ^^^C is transferred with respect to 
^^C (the predicted isotopic ratio can be as high 
as 10'*, hence lO'' ^•^C nuclei are transferred for 
a single ^''C nucleus). The amount of ^^^C in 



^* Note that the values of C, N, Na, C isotopic ratios in the 
following discussion are those listed in Table [2] and |4] 
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the atmosphere of the CEMP star before it be- 
comes CEMP counts for small amounts of mass 
transfer. Fig. [6] shows the predictions for the 
same models shown in Fig. \5\ under the as- 
sumption that the original (pre mass-transfer) C 
abundance and C isotopic ratios in the observed 
stars were [C/Fe]=-0.4 and 20, which is a typ- 
ical value measure d for m etal-poor giants (see 
e.g. iGratton et al.l (|2000n ). The models in the 
plot take into account different degrees of dilu- 
tion, corresponding to different positions along 
the lines. The C isotopic ratios measured in 
our sample are all quite low, some almost at the 
equilibrium value, much lower than predicted by 
low-mass shell-nucleosynthesis models (^10,000 
see e.g. Cristallo et al. 2009), while the ratios 
for more massive stars (>2.5Mq) are closer to 
the measured ones. However, a star polluted by 
a star with a mass >2.5M0, where HBB oper- 
ates should be characterized by [N/C]>0, which 
is only found in two of the stars in our sample, 
making this possibility quite unlikely (see also the 
discussion in the next section). 

During their evolution, the stars in our sam- 
ple undergo several mixing episodes (i.e. first 
dredge-up, extramixing corresponding to the 
RGB bump) from the main sequence to the evo- 
lutionary state we observe them in today. This 
mixing episodes are known to affect the C iso- 
topic ratio in C normal stars. It is reasonable to 
assume that they also might affect the accreted 
C, altering the C isotopic ratio of the accreted 
material. 

This effect is however likely not important. 
In fact, similar values for C isotopic ratios 
are found in uneyolved CEMP-s stars (see 
e.g. iLucatello et al"]|2003l : iMasseron et al.|[2010al : 
lAoki et al.l 120021 ). which have not undergone any 
of such mixing episodes and have thus maintained 
the original ^^C/^'^C in transferred material. 
Moreover, the extent of the effect of the mixing 
episodes on the accreted C isotopic ratio should 
likely depend on the time in its evolution when 
the C-rich material is transferred. Since it is 
sensible to assume that the CEMP-s stars in 
our sample underwent mass transfer at different 
times in their evolution, the fact that all the mea- 
sured isotopic ratios are all quite low, argues that 
the effect of the mixing episodes on such quantity 
should be marginal. Therefore, it is reasonable to 
conclude that the measured C isotopic ratios in- 
deed reflect that in the polluting material. 

To account for such a high N abundance and low 
C isotopic ratios in low-mass AGB stars, we con- 
firm that some degree of of extra-mixing with a 
partial activation of proton capture nucleosynthe- 
sis, the so-called Cool Bottom Processin g (CBP), 
needs to be invoked (jNollett et al.l 12003 ). How- 
ever, in most published works, such process is 
treated as a parameter and thus the amount of 
CN processing is highly uncertain, making the 
individual C and N abundance predictions much 
more sensitive to modeling details than their sum. 
This process is known to occur during the first red 



giant phase in stars with M<2.5 M© ([Charbonnel I 
Il995n and was sugge sted to occur between dredge 
up in AGB stars ([Nollett et al.l l2003h . It con- 
sists of circulation of material from the base of 
the convective envelope into the thin radiative re- 
gion located on top of the H-burning shell, where 
the material is processed by proton captures and 
then brought back to the envelope, producing 
the signature of CN processing at the stellar sur- 
face. The low 12C/13C ratios, high N (but with 
[C/N]>0) observed in CEMP-s stars indicates 
that CBP also takes place in low' metallicity AGB 
stars (see also lAoki et al.l 120021 ; IMasseron et ahl 
l2010al ). Unfortunately, the physical mechanism 
responsible for such process is still not clear. Sev- 
eral scenarios have been proposed, e.g., the mix- 
ing process could be driven by a molecular weight 
invers ion created by^he ''He('3He,2p)'*He reaction 
(Egg leton et all [20081 ). or by magnetic buoyancy 
at the bottom of the envelope (e.g., Palmerini, 
Nollett, &: Busso 2008). 

While we obtained measurements for O for all 
the stars in our sample, a comparison with the 
model predictions would not be meaningful. All 
the models discussed in the text start with solar- 
scaled composition. While this is not impor- 
tant for the other elements we discuss as their 
production in the AGB nucleosynthesis is much 
larger than the initial content, this is not true 
for O, which is in fact typically produced at the 
[O/Fe]'^0.5 level, a value comparable to that mea- 
sured in metal-poor stars, making the adoption of 
a correct initial O for the models important. The 
typical predicted values for O are in fact lower 
than the measured ones; however it is presently 
not possible to determine whether this is due to 
the incorrect assumption on the model star ini- 
tial composition or issues in the understanding of 
O nucleosynthesis in low-metallicity AGB stars 
(see for instance the impact of overshooting on O 
abundance, e.g., Herwig 2005 ARA&A). 

5.3. Comparison with observations: F 

F was measured only in two CEMP stars in our sam- 
ple, both s-process rich, while for the remaining eight 
objects (including the two CEMP-no stars) only upper 
limits could be placed. Therefore, as far as ^^F is con- 
cerned, no meaningful comparison is possible between 
CEMP-s and CEMP-no stars, and not much informa- 
tion on the formation scenario of CEMP-no stars can be 
derived from the present data. 

Fig. [7] shows the derived F abundances with respect to 
the sum of the C and N abundances in the sample stars. 

Results from four different families of models are plot- 
ted for comparison to the data; they provide yields for 
of various elements among which are C, N, Na, and F. 
In Fig. [3 for each individual model different degrees of 
dilution correspond to different positions along its corre- 
spondent line. 

In the upper panel the model plotted are all com- 
puted for a mass of 2M0 and a metallicity of 
[Fe/H] — — 2.3. The dotted line indi cates the predic- 
tion from Karakas fc Lattanziol (|2007f) . the one used in 
Xugaro et al.. (,200Sl to account for the F content re- 
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ported in iSchuler et al.l (I2008D . The solid line describes 
the model from iKarakas et all (|2008f) . which includes a 
^^C pocket of O.OO2M0 at the top layer of the He- and 
C-rich inter-shell region and adopts the upper limit for 
the rate of the ^^F(a, p)^^Ne reaction (in an attemp to 
reproduce the high F abundances in the Jorissen sam- 
ple). Among the models available in the literature, this 
is the one that for the mass and metallicity under con- 
sideration, produces the largest amount of F with re- 
spect to C and N. The dot-das hed line indicates the pre- 
diction of the model from Cristallo et alj (|2009f ).. which 
accounts for the appropriate C-bearing molecules opaci- 
ties in the stellar evolution model, unlike the two other 
models, which adopt solar scaled compositio n. Finally, 
the dotted line shows the prediction from the iStanclifFel 
(EOOa) model, which uses solar-scaled opacities. 

As shown in the upper panel of Fig. [71 the 
two F measurements can be accounted for by the 
IKarakas fc Lattanziol ()2007[ ) or the Cristallo et al. (2009) 
2 M0 star models. On the other hand, the F predicted 
from the Karakas et al. (2008) 2 M0 star is higher. All 
the upper limits excepting one derived for F in our sam- 
ple fall lower than the two measured abundances. Two 
of the upper limits are in agreement with the predictions 
of Stancliffe & Glebbeek (2008) in the figure. We note 
that a significant fraction of the upper limits are located 
under the predicted lines. 

As mentioned before, for the Lugaro et al., Cristallo et 
al. and Karakas et al., F production peaks at ^ 2.2 Mq 
(see iLugaro et al.l I2008D . Quantitative predictions are 
provid ed o nly for 2.0 M(7) A GB models in lCristallo et al.l 
(pool and IKarakas et al.l (pOOSi ) papers, so no compari- 
son with different mass models could be perform ed. On 
the other hand, iStancliffe &: Glebbeekl (2008) predic- 
tion is lower than the measurements, possibly account- 
ing for two of the upper limits (one CEMP-s and one 
CEMP-no star). However, it still cannot account for 
most of the upper limits derive d for CEMP-s stars. In the 
IStancliffe fc Glebbeek I (|2008[ ) models, unlike the others, 
F production peaks at ~1.5Mq. Hence, we cannot in- 
voke lower mass models to reconcile the discrepancy in 
t his case. 

IKarakas fc Lattanziol ()2007l ) provide F production pre- 
dictions for several AGB masses. In the lower panel 
of Fig . [7] the predictions for the IKarakas fc Lattanziol 
([20077 niodels are plotted together with the data. Mod- 
els are for masses 1.25, 1.75, 2.0, 2.5, 3.0, 3.5 and 4Mq 
at a metallicity of [Fe/H]= —2.3. The models at 1.25, 
1.75 and 2.5 Mq seem to reproduce satisfactorily the two 
F measurements in our sample, even though it should 
be kept in mind that, as mentioned, their metallicity is 
higher than that of the models ([Fe/H]=— 1.3 and —2.0). 
However, as discussed later in this section, the effect of 
metallicity on [F/C-l-N] is likely very small. On the other 
hand, when considering the measurements obtained for 
the whole sample, it can be noticed that predictions for 
F to C-l-N for masses from 1.25 to 2.5 Mq are higher than 
four out of the eight CEMP-s stars. While the predic- 
tions for models for 3, 3.5 and 4 Mq can account for the 
upper limits F to C+N ratios in those four CEMP-s, it 
is important to keep in mind that for [Fe/H] < —2 all 
AGB stars with intial mass higher than M~2.5Mq will 
be affected by HBB, which burns C producing N (e.g.. 



Karakas fc Lattanzio 2007). A star polluted by such an 
object, would not be a C-rich, but rather N-rich, which 
is not the case of our sample (three out of those five 
CEMP-s have [C/N]>0). Moreover, stars of this mass 
range are not very efficient in producing s-process ele- 
ments (which are enhanced in eight stars in our sample, 
we provide Ba abundances in Table [5] for completeness). 
Finally, a simple back-of-the-envelope calculation adopt- 
ing a Salpeter IMF and assuming that, in the considered 
mass range, the likelihood of having a companion is not 
affected by mass, indicates that the probability of observ- 
ing, among a sample of eight CEMP-s stars with a low- 
mass companion (1.25<M<4M0), four object that had 
a companion in the 2.5 to 4M0 mass range, is very small 
(lO"'^). For these reasons, the possibility that the abun- 
dance ratios observed are due to stars more massive than 
2.5 Mq seems highly unlikely. We note that for one 
of the CEMP-no stars, CS 30314-67, the derived 
upper limit to F abundance is also quite low and, 
similarly to most CEMP-s stars in our sample, is 
not reproduced by the low-mass, low-metallicity 
AGB models. This could be in principle used 
as an argument against the scenario which at- 
tributes their peculiar composition to pollution 
from low-mass AGB. However, given the issues 
that the models have in reproducing the obser- 
vations as far as F is concerned in CEMP-s stars 
(which are known to arise from AGB mass trans- 
fer), the AGB pollution scenario cannot be ruled 
out. 

Notice that for each model discussed we consider the 
last thermal pulse. According to e.g. Izzard fc Tout 
(2003) AGB evolution can be truncated by binary in- 
teraction, with the AGB star experiencing fewer ther- 
mal pulses with respect to predictions of an isolated ob- 
ject. To illustrate the evolution of the abundance ra- 
tios through different thermal pulses we plot log e(F) and 
loge(C-hN) from 2Mq (Cristallo et al. 2009), overplot- 
ting our measurements for comparison (see upper panel 
of Fig. [8]). With the exception of one limit, in principle 
most of the upper limits could be explained by predic- 
tions of this particular model at different thermal pulses. 
In the lower panel, on the other hand, we show predic- 
tions from Cristallo et al. (2009) for a 2 Mq with two dif- 
ferent metallicities, [Fe/H]=— 1.3 and —2.3. The figure 
shows a change of an order of magnitude in metalicity 
has only a small effect on the ratio [F/C-l-N]. For this 
reason, even though our sample star cover a large range 
in metallicity, we compare our measurement to the pre- 
dictions at a single metallicity, the lowest one available. 
We invite the reader to keep in mind that this is an ex- 
trapolation to lower metallicity of the behaviour of the 
model and as such quite uncertain. 

Beside the sources of uncertainties discussed in the pre- 
vious paragraphs, as mentioned before, currently AGB 
nucleosynthesis calculations do not include GBP. Indeed, 
at typical AGB envelope conditions F can be depleted 
significantly by extra-mixing processes, depending on 
their efficiency (Palmerini, NoUett, fc Busso 2008), via 
the ^^F(p,q;) reaction. Therefore, the inclusion of an 
appropriate treatment for CBP in AGB modeling could 
possibly decrease the F production, bringing the predic- 
tions closer to the observations. It is noteworthy that 
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partial F depletion has been proposed to be active at 
the tip of the previous Red Giant Branch (RGB) phase 
(Denissenkov et al. 2006). It must be kept in mind, 
however, that such a mechanism has been invoked to ex- 
plain the decrease of F abundance in stars at the tip of 
the RGB in the globular cluster M4, which, as likely 
all globular clusters, show abundance v ariations which 
are n ot observed in field stars (see e.g. iCarretta et al.l 
|2009() . It is also important to note that one of the two 
stars in our sample for which we can actually measure F, 
HE 1152-0355, is quite close to the RGB tip. 

5.4. Comparison with observations: Na 

Models for low-mass AGB stars predict quite success- 
fully the observed Na abundances. Na is produced by 
AGB stars in different ways. First, during the TP 
neutrons produced via the ^^Ne(Q!,n)^^Mg may be cap- 
tured by the ^^Ne itself, producing ^'^Na via the reaction 
chain ^^Ne(n,7)^^Ne(/3~)^^Na. Secondly, sodium may 
be strongly produced also by ^^Ne(p,7)^'^Na. Indeed, 
after the TDU, at similar mass coordinate of the ^''N- 
pocket a ^^Na-pocket is formed in radiative conditions, 
where ^'^Na is given by the equilibrium between the pro- 
duction channel ^^Ne(p,7)^^Na and the depletion chan- 
nels ^3Na(p,7)^''Mg and 23Na(p,a)2°Ne (see Cristallo et 
al. 2006 for details, and Bisterzo et al. 2006). Finally, 
a marginal contribution is given by neutron capture on 
^^Ne in the ^^C-pocket. Despite these different processes, 
Na production in AGB stars is less affected by nuclear 
uncertainties compared to F. For these reasons, AGB pre- 
dictions for Na are expected to be more robust than F. 
Fig. |9] shows Na and C-l-N in the sample stars together 
with the predictions from the same models as in Fig. [71 
The upper panel shows two of the models represented in 
the correspondent panel of Fig. [7] because those are the 
only ones that provide Na abundance predictions. It is 
noteworthy how the discrepancy in Na between these two 
models (~0.05 dex at a given log e (C-l-N) is much smaller 
than that in F between the same models (~0.5 dex), sup- 
porting the idea that Na production is generally bet- 
ter understood and/or less sensitive to the recipes and 
parameters adopted in different families of models. As 
can be noticed from in Fig. 1101 the measured Na abun- 
dances are consistent with pollution arising from a 1.25 
to 1.75 Mq AGB star. Six out of seven sample CEMP-s 
for which Na was measured and both CEMP-no stars are 
consistent with having a companion in that mass range. 

6. CONCLUSIONS 

We have obtained IR high-resolution observations in 
a sample of ten CEMP stars, eight CEMP-s and two 
CEMP-no. The aim was to measure Na abundances, ^^F 
contents from HE, O from CO and N from CN lines in 
the K band. 

A comparison of predictions from different families of 
low-metallicity, low-mass (1.25 to 1.75 M© range) AGB 
nucleosynthesis models show that they reproduce well 
the observed Na abundances in the CEMP-s stars in our 
sample. This result strongly argues in favor of a polluter 
of low-mass for these objects, in agreement with the cur- 
rently most accepted scenario for their formation. 

F could only be measured in two CEMP-s stars, while 
for the remaining eight objects only upper limits could be 



derived. However, these data are sufficient to show that a 
range of [F/C-l-N] values are produced in low-metallicity 
AGB stars, in accord with predictions of the mass- and 
metallicity-dependence of F production in AGB models. 

The only two F measurements obtained (the CEMP- 
s stars HE 1152-0355 and HD5223 respectively at 
[Fe/H]=— 1.27 and —2.06) are accounted for, within the 
errors, low-mass, low-metallicity AGB models. On the 
other hand, most of the derived upper limits for F abun- 
dance in CEMP-s are not satisfactorily accounted for 
by nucleosynthetic computations. In fact, the compar- 
ison with four of the most recent models for low-mass 
(2.0 Mq), low-metallicity AGB nucleosynthetic models 
shows that there are large differences in the predictions 
between different families of models, which cannot re- 
produce several of the upper limits, not providing any 
[F/C-fN] ratios predictions low enough to account for 
the values meas ured in several of the samp le stars. A 
comparison with lKarakas fc Lattanziol ()2007f) models for 
different masses indicates that only objects more mas- 
sive (2.5 M0 <M<4M0) that those generally consid- 
ered as responsible for CEMP-s peculiar abundances 
(M~1.5M0) produce the [F/C-l-N] ratios observed in 
about two thirds of our CEMP-s sample. This possi- 
bility is not only unlikely because of simple IMF con- 
siderations, but it is also challenged by the fact that 
at low metallicity, HBB is active in stars in this mass 
range, producing N and depleting C (b ringing the ratio to 
[N/C]^l (see e.g. iJohnson et ahl 12007 )) and preventing 
the star from becoming a CEMP-s polluter. As shown 
in FiglH, the lKarakas fe Latta nzio (2007) 1.25 Mq 
model is the low-mass AGB ones that comes clos- 
est to explaining the upper limits. Given the 
discussed spread present in F production in dif- 
ferent models, it is indeed possible that it af- 
fects also predictions for other masses, including 
1.25 Mq , hence accounting for the derived up- 
per limits. It is however impossible, at present 
time, to systematically test this differences, as 
only Karakas & Lattanzio (2007) provide F pro- 
duction calculations for several AGB masses. 

We discussed two possible solutions that could explain 
the lowest F upper limit range. One is that the standard 
evolution of an AGB star may be truncated by binary 
interaction, and a lower F abundance can be obtained in 
early TP pulses (even though this would lead to lower 
Ba production, which could make the accounting for ob- 
served Ba abundances possibly more challenging). The 
other is the action of CBP, which may reduce the F in the 
AGB envelope if such extra- mixing processes exposes ma- 
terial at temperatures high enough to activate ^^F(p,a). 

The conclusion that can be drawn from the compari- 
son between the data and the current models is that F 
production in low metallicity AGB stars is probably not 
as high as expected on the basis of the current models. 
It is important to keep in mind that, as discussed in the 
introduction, at solar metallicity model predictions are 
in agreement with the observations, therefore such prob- 
lem seem to be peculiar to the low metallicity regime. 
AGB nucleosynthesis, especially at low metallicity, is 
still not fully understood and a large number of uncer- 
tainties affect in particular F nucleosynthesis. Improve- 
ments in nuclear reaction rate (such as ^'^C (a, 7)^*0 and 
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^^F(a,p)^^Ne) accuracy are needed, and a better grasp 
of the CBP mechanism is highly desirable. O nly one of 
the available models, bv lCristallo et all (|2009[ ) takes into 
account the appropriate C, N and O enhancements in the 
computation of the stellar opacities; also, the prescrip- 
tion for mass loss is still rather uncertain. Moreover, 
currently no model takes into account the effect on the 
AGB evolution of the presence of a (close) companion, 
which could stimulate extra mixing and/or truncate the 
AGB phase (see e.g. Izzard & Tout 2003), affecting the 
nucleosynthesis. The inclusion of these ingredients into 
a new generation of AGB models might in the future be 
able to account for the present observations. 

From an observational point of view, pursuing a larger 
number of F abundance measurements in CEMP stars 
is extremely important. In particular, K-band high- 



resolution spectral observations of cool (Te// <4200K) 
CEMP stars would allow actual measurements (rather 
than upper limit placements) of F even for low F abun- 
dances (loge(F)~2). These measurements would provide 
more stringent constraints that are urgently needed for a 
future generation of AGB models and to effectively probe 
F chemical evolution. 
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the National Research Council of Canada. 
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Fig. 1.— HF line at 2335.8 nm in star HD 5223 with synthetic spectra for log e(F)=3.95, 4.0 and 4.05. 
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Fig. 2. — Measured F abundance as a function of T^g in the sample stars. Inverted triangles indicate upper limits. There is a clear 
correlation between the upper limits set on the F abundances and the Teg , indicating that the most stringent limits can be obtained from 
cooler stars. 
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Fig. 3. — Comparison of the derived abundances and literature for N, O and Na. Sources for literature data as in Table [2] 
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Fig. 4. — Example of comparison between a spherical model (solid black line) and a plane-parallel one (red dashed line). The models are 
computed for the atmospheric parameters of HE 1152-0355 
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6.5 7 7.5 8 8.5 

log e(C + N) 



Fig. 5. — Upper panel: N abundance as a function of C+N in 10 CE MP stars. Green symb o ls are CEMP-no stars. Theoretical modeling 
results are indicated by lines. Upper panel; plotted models are from iKarakas &: Lattanz io (2007) for different masses. Solid black line 
1.25 M0, red dotted line 1.75 M©, teal dashed line 2.25 Mq, green dot-dashed 2.5 Mq, magenta dot-dashed line 3.0 Mq and short-long 
dashed blue line 4Mq. Lower panel: Model predictions from Cristallo et al. (2009) for a 2Mq star at a metallicity of -1.3 (dotted teal 
line) and -2.3 (dashed gold line). 
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Fig. 6. — Log arithm of C isotopic r a tio pl o tted as a function of C -l-N. Upper panel: the oreti cal modeling resu lts are indicated by lines. 
Prediction from IKarakas fc Lattanziol I I2007I ). iKarakas et all II2008I ). ICristallo et ahl II2009I) and 'Sta nch^ fcoog ') are plotted. The models 
are indistinguishable in the plot as they all predict very similar (high) C isotopic ratios. All models are for M=2M0 and [Fe/H]=— 2.3. 
Symbols in the lower panel are as in Fig. \5\ upper panel. Note that in this case the initial abundance (pre mass transfer) of the observed 
stars for -'^'^C is taken into account (see text for details). 
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Fig. 7. — Upper panel: F abundance as a function of C+N in 10 CEMP stars. The inverted triangles indicate upper limits, green symbols 



are CEMP-no stars. Theoretical modeling results are indicated by lines. Dotted line indicates the prediction from 'Karakas & Lattanzic 
2007), the solid line is from Karakas ct al. (2008), the dot-dashed line is from[Cristallo ot al. (2009) and the dashed line is from Stancliff^ 
2009). All models are for M=2Mq and [Fe/H]=— 2.3. Lower panel: as in upper panel, but plotted models are from lKarakas &: Lattanzic 



(2007) for different masses. Solid black line 1.25 Mq, red dotted line 1.75 Mq, teal dashed line 2.25 Mq, green dot-dashed 2.5 Mq, magenta 
dot-dashed line 3.0 Mq, short-long dashed yellow line 3.5 Mq and short-long dashed blue line 4Mq. 
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Fig. 8. — Symbols for measurements in the sample are as in Fig. [7] Upper panel: log e(F) and log e(C-l-N) for a 2 Mq star at [Fe/H]=-2.1 
(Cristallo et al. 2009) for thermal pulses 1 (black solid line), 3 (dotted green line), 5 (short-dashed red line), 10 (long-dashed blue line) 
and 15 (last pulse, dot-dashed magenta line). Lower panel: Model predictions from Cristallo et al. (2009) for a 2Mq star at a metallicity 
of -1.3 (dotted teal line) and -2.3 (dashed gold line). 
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Fig. 9. — Na plotted as a function of C+N. Symbols as in Fig. [7] Note that only two of the four models plotted in the upper panel of 
Figure [7] are plotted h ere as they are the ones reporting predictions for Na production. The lower panel shows Na production prediction in 
ICristallo et"all l(2009h at [Fe/H]=-1.3 and -2.3. 
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TABLE 1 
Observational log 



Star 




RA 






Dec 


V 


K 


Obs 




Exp. 


ID 




















Date 




Time (s) 


CS 22942-19 


00 


57 


18 


136 


-25 


26 


09.69 


12.69 
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2007- 


10- 


21 


1200 


CS 29497-34 


00 


41 


39 
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-26 


18 


54.48 


12.22 


9.996 


2007- 


10- 


26 


900 


CS 30314-67 


20 


52 


50 


986 


-34 


19 


40.52 


11.85 


9.028 


2007- 


10- 


23 


900 


CS 22948-27 


21 


37 


45 


774 


-39 


27 


22.31 


12.660 


10.427 


2007- 


10- 


27 


1200 


CS 29502-92 


22 


22 


35 


995 


-01 


38 


27.53 


11.87 


9.601 


2007- 


10- 


24 


1200 


CS 30322-23 


21 


30 


32 


184 


-46 


16 


24.72 


12.21 


9.313 


2007- 


10- 


23 


900 


HE 1152-0355 


11 


55 


06 


057 


-04 


12 


24.71 


11.43 


8.429 


2008- 


02- 


26 


720 


HE 1305-1-0007 


13 


08 


03 


849 


-00 


08 


47.40 


12.223 


9.600 


2008- 


02- 


27 


1200 


HD 187861 


19 


56 


26 


945 


-65 


22 


08.01 


9.19 


6.699 


2007- 


10- 


21 


180 


HD 5223 


00 


54 


13 


610 


-f24 


04 


01.51 


8.47 


5.673 


2007- 


10- 


21 


180 


HD 122563 


14 


02 


31 


846 


-f09 


41 


09.94 


6.20 


3.731 


2008- 


02- 


28 


45 



TABLE 2 

Adopted atmospheric parameters, C and Ba abundances 



Star 


TefF 

(K) 


logs 


[Fe/H] 


[C/Fe] 


[Ba/Fe] 


ID 




(dex) 


(dex) 


(dex) 



CS 22942-19 


5100^1100° 


2.5±0.1° 


-2.5^10.15" 


2.2±0.10'" 


1.76±0.10'' 


CS 29497-34 


4800^= 


1.8^= 


-2.9^= 


2.6'= 


2.03'= 


CS 30314-67 


4400±100'' 


0.7±0.3'' 


-2.85±0.18'' 


0.8±0.18'''' 


-0.57±0.14'' 


CS 22948-27 


4800^= 


1.8^= 


-2.5^= 


2.4= 


2.26'= 


CS 29502-92 


4890±100'* 


1.7±0.2'' 


-3.18±0.09'* 


0.96±0.22'''» 


-1.26±0.09<* 


CS 30322-23 


4100±100'= 


-0.3±0.3'= 


-3.39±0.18'= 


0.6"= 


0.54±0.10'= 


HE 1152-0355 


4000±90-^ 


1.0±0.3-^ 


-1.27±0.27^' 


0.58±0.2-'' 


1.58-'' 


HE 1305-1-0007 


4560±909 


1.0±0.39 


-2.59 


2.4±0.359 


2.9±0.59 


HD 187861 


4600±150" 


1.7±0.2" 


-2.36^10.21" 


2.02±0.10'" 


1.39 ±0.10" 


HD5223 


4500±90-^ 


1.0±0.3-^ 


-2.06±0.13^ 


1.57^ 


1.82^ 


HD 122563 


4600±100'' 


1.1±0.3'' 


-2.82±0.10'' 


-0.10±0.05'''' 


-1.28±0.23 



Not e. — alMasseron et al.l j2010bl). blAoki et al.l ll2002f) . clBarbuv et al.N2005ft . dILai et al. 
20071) . e lMasseron et al.l I I2006I) . f lGoswami et al.l 120061) . g lBeers et al.l I I2007I) h ISpite et al. 
200|), i C abundance has been adjusted (see text for details) 



TABLE 3 

Sensitivity of the derived abundances to the atmospheric parameters for 
HE 1152-0355 and HD 5223 



Element 


-(-200 K 


HE 1152-0355 
logg [Fe/H] 
-1-0.5 dex -1-0.3 dex 


[C/Fe] 
0.3 dex 


-1-200 K 


HD 5223 
log 3 [Fe/H] 
+0.5 dex -1-0.3 dex 


[C/Fe] 
0.3 dex 


O 


0.4 


0.1 


0.1 


-0.1 


0.5 


0.1 


0.1 


-0.1 


N 


0.4 


0.1 


0.1 


-0.2 


0.4 


0.1 


0.1 


-0.2 


F 


0.5 


0.0 


0.2 


0.1 


0.5 


0.0 


0.2 


0.1 


Na 










0.2 


0.1 


0.1 


0.1 
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TABLE 4 

N,0, Na and F abundance measurements for, the program stars. C isotopic ratio is 

ALSO listed 



Star 
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Star 




log 6(F) 


[F/Fe] 


a(loge(F)) 


loge(Na) 


[Na/Fe] 


f7(loge(Na)) 


CS 22942-19 
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Note. — The adopted solar abundances arc from|G revesse et all l2007l) for C,N, O, F and Na, and 
IGrevesse fc Sauvaj ilQQgT for Ba. 
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